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摘 要 


高 负荷 跨 音 透 平 导 叶 中 ,流动 过 程 十 分 复杂 ,DDES 方法 在 该 条 件 下 对 流动 细节 地 捕 提 具有 特殊 优势 . 本 文采 


用 自主 开发 的 DDES 程序 , 获得 可 靠 而 精确 的 计算 结果 , 研究 跨 音 叶片 内 部 的 流动 机 理 . 根据 不 同 区 域 的 流动 现象 , 作者 
分 别 研究 叶片 吸力 面 处 , 激 波 和 边界 层 的 相互 作用 ; 叶片 尾 流 区 , 激 波 和 尾 迹 的 相互 作用 . 并 对 这 些 区 域 进 行 粹 生成 情况 


分 析 , 为 后 续 跨 音 透 平 中 的 损失 减 小 提供 理论 基础 . 


关键 词 DDES; 跨 音 透 平 导 叶 ; 激 波 ; 边界 层 ; 尾 迹 
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Key words 


0 引 言 


燃气 轮机 负荷 正 朝 着 越 来 越 高 的 方向 发 展 ， 脱 
胀 比 逐 渐 增 大 , 在 透 平 级 数 不 变 的 条 件 下 , 这 势必 会 . 
导致 单 级 负荷 的 增加 ， 甚 至 会 在 通道 内 出 现 跨 

音 流 动 .通道 中 激 波 诱导 的 边界 层 分 离 以 及 诗 片 
尾部 的 脱落 涡 等 一 直 是 影响 效率 提高 的 主要 因素 ; 
采用 数值 模拟 的 方式 , 可 以 在 降低 实验 成 本 的 同时 ， 
获取 更 多 的 流 场 信息 ， 这 对 人 燃 机 内 部 流动 的 物理 分 
析 十 分 重要 ，. 

对 于 燃 机 中 复杂 的 清流 流动 , LES 能 够 得 到 比 
RANS 更 精细 的 模拟 结果 ,而 RANS 消耗 的 计算 资 
源 远 小 于 LES. 为 此 ,1997 年 Spalart 等 中 提出 了 
将 RANS 和 LES 混合 的 DES(Detached Eddy Sim- 
ulation) 方法 , 在 保证 精度 的 同时 ,降低 对 计算 资源 
的 需求 。2006 年 ，Spalart 等 争 修改 交界 处 混合 方 
法 的 判定 准则 , 在 近 壁 区 延迟 LES 模拟 , 被 称 为 延迟 


收 稿 日 期 :2016-12-27; 修订 日 期 : 2017-04-20 


DDES; transonic turbine guide vane; shock wave; 


e and wake at wake region of the blade. And the entropy generation 
hich provides the theoretical basis of the los ‘duction in transonic turbine. 


undary layer; wake 


© 


\ 
有 Sy DDES(Delayed DES), 修正 了 DES97 中 
y 为 模 化 不 足 的 问题 , 使 DDES 能 够 和 实验 更 好 地 
吻合 。 因 此 ， 在 激 波 处 ， 边 界 层 分 离 区 ， 尾 迹 涡 脱 
落 等 复杂 流动 区 域 采 用 DDES 可 以 在 获得 精细 流 
场 的 条 件 下 ， 减 小 计算 量 ， 适 合 跨 音 导 叶 内 部 机 理 
研究 . 

本 文 使 用 完全 自主 开发 的 多 块 并 行 求解 器 加 ， 
采用 多 块 并 行 计算 以 增 大 计算 规模 ,利用 双 时 间 步 
法 以 达到 透 平 流动 中 非 定常 求解 要 求 ; 对 激 波 ,， 边 
界 层 和 尾 迹 等 复杂 流动 区 域 ， 程 序 采 用 耗 散 项 自 适 
应 函数 四 , 提高 计算 精度 ,获得 较 好 结果 . 


1 研究 对 象 及 方法 
1.1 求解 域 及 网 格 
本 文 计算 几何 模型 选取 的 是 加 拿 大 国家 研究 委 
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员 会 (CNRC) 测试 使 用 的 高 负荷 透 平 进口 导 叶 回 ， 
计算 工 况 为 跨 音 工 况 ， 可 用 于 研究 激 波 和 边界 层 以 
及 尾 迹 的 相互 作用 , 为 除去 上 下 端 壁 对 流动 的 干扰 ， 
选取 导 叶 中 截面 几何 进行 模拟 , 为 进一步 简化 , 选择 
平面 叶 栅 的 流动 状况 进行 讨论 ， 叶 型 表面 设置 为 绝 
热 壁 面 . 主要 数据 由 表 1 给 出 . 选取 的 计算 域 由 图 1 
给 出 ,入口 距 叶片 前 缘 选 取 0.72 个 弦 长 ， 出口 距 叶 
片尾 缘 选 取 1.5 个 弱 长 , 周 向 长 度 选取 一 个 栅 距 , 展 
向 计算 区 域 高 度 选 定 10 mm， 沿 展 向 和 周 向 均 设 置 
为 周期 性 条 件 . 


图 2 尾 缘 网 格 
Fig. 2 Trailing edge mesh 


表 1 叶 型 参数 
Table 1 The parameters of the blade 


弦 长 /mm 轴 向 弦 长 /mm ” 栅 距 /mm 展 长 /mm 9 
204 84.7 147.8 112.8 


结果 分 析 
2.1 ,说 扯 结 果 与 实验 数据 对 比 
Ce 3.5 个 通过 周期 后 ， 取 叶片 尾 缘 监 测 点 压 
~、 邹 时 间 的 变化 做 信里 叶 变 换 ， 做 功率 谱 密 度 分 析 
AN 为 图 3， 峰值 处 频率 为 25.92 kHz， 满 足 文献 实验 给 
定 日 的 涡 脱落 频率 参考 范围 24~26.5 kHz。 
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图 1 计算 域 


Fig. 1 Computation domain 
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将 叶 型 几何 数据 导入 NUMECA 中 ,利用 Auto \\ \ 1 
Grid5 生成 满足 要 求 的 计算 区 域 ， 再 导入 IGG 中 . 十 从 局 人 

行 分 类 ,共计 95 块 . 为 保证 计算 的 精度 ,全 


:hin 
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结构 化 网 格 ， 且 在 网 格 数 设置 时 ， 保 证 计算 过 程 中 
可 以 采用 多 重 网 格 方法 进行 加 速 收敛 ; 加 密 后 ， 网 
格 总 数 为 1.4x10"” . 叶片 尾 缘 采 用 圆 弧 拟 合 , 保证 尾 
缘 区 域 的 光滑 性 ， 尾 缘 部 分 网 格 如 图 2 所 示 . 在 计 
算 中 , 时 间 步 长 选取 10  s, 每 个 物理 时 间 步 欠 代 结 
束 的 条 件 为 内 达 代 30 步 或 残 差 下 降 4 个 量 级 . 
1.2 边界 条 件 

音 计算 工 况 , 出 口 马 赫 数 选取 1.16, 其 他 进出 
口 参数 如 表 2 所 示 : 


表 2 进出 口 参数 
Table 2 Inlet and outlet parameters 


进口 总 温 /K 进口 总 压 /Pa 膨胀 比 。 来 流 攻 角 /(°) 
293 101300 2.3 —10 


Frequency 


图 3 功率 谱 密度 图 


Fig. 3 Power spectral density 


2.2 能 激 疲 分 析 

做 出 瞬时 密度 梯度 图 (DGM, Density Gradient 
Magnitude), 采用 流 场 中 最 大 密度 梯度 进行 单位 化 ， 
得 图 4, 其 中 叶片 尾 缘 , 两 道 主 斜 激 波 之 前 , 还 存在 
两 道 唇 部 激 波 ， 局 部 放大 后 如 图 5 中 左 图 所 示 ， 层 
部 激 波 强度 较 弱 ,产生 原因 为 尾 缘 圆 弧 直 径 稍 大 于 
尾部 厚度 ， 

此 外 ,图 4 瞬时 密度 梯度 图 中 ， 在 吸力 面前 段 
部 分 , 出 现 一 个 密度 梯度 较 高 的 区 域 , 局 部 放大 后 如 
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图 5 中 右 图 所 示 ， 对 时 均 结 果 做 马赫 数 云图 得 图 4 
右 图 , 黑 线 是 马赫 数 为 1 的 等 值 线 , 所 以 图 5 中 , 在 
吸力 面前 段 密 度 高 区 为 正 激 波 ， 产生 原因 为 本 次 计 
算 选 用 叶 型 内 部 流 道 为 渐 缩 流 道 ,吸力 面 膨 胀 较 快 ， 
当 流速 达到 音速 时 ， 后 续 流 道 无 法 为 气动 提供 必要 
的 膨胀 条 件 , 所 以 出 现 正 激 波 , 在 叶片 压力 面 处 , 因 
为 流体 膨胀 较 慢 ， 故 相同 流 道 截面 位 置 ， 并 未 出 现 
该 现象 . 


瞬时 密度 梯度 时 均 马 ag 
图 4 弱 激 波 


Fig. 4 Weak shock ee 


和 
7 


图 5 弱 激 波 局 部 放大 图 


Fig. 5 Partial enlarged view of weak shock wave 


正 激 波 


2.3 激 波 与 边界 层 干涉 机 理 分 析 

在 激 波 边界 层 作 用 区 域 ， 对 时 均 马 赫 数 云图 进 
行 局 部 放大 得 图 6, 可 以 看 到 斜 激 波 在 叶片 边界 层 发 
生 反 射 , 同时 诱导 边界 层 分 离 。 做 出 该 区 域 的 形状 因 
子 , 结果 表明 , 在 吸力 面 处 出 现 边 界 层 分 离 现 象 , 并 
产生 分 离 泡 ， 用 距 壁 面 的 相对 位 置 表征 干涉 区 形状 
因子 的 大 小 ,可 以 看 到 , 干涉 区 形状 因子 突然 增加 ， 
也 说 明 该 区 域 发 生 了 边界 层 分 离 。 


NN 
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: \ N | Mach 


图 6 边界 层 分 离 


Fig. 6 Boundary layer separation 


Tb 图 7(a) 
为 湾 方 能 图 , 在 两 道 激 波 处 , 激 波 交 汇 处 以 及 分 离 泡 


已 


人 速度 波动 较 大 . 做 出 马赫 数 云图 ， 
壁 区 用 黑 实 线 表示 马赫 数 为 1 的 等 值 线 ， 壁 面 附 


近 黑 色 区 域 为 诱导 产生 的 分 离 泡 ， 用 线 图 简化 干涉 
过 程 . 激 波 前 后 速度 波动 是 该 处 注 动 能 较 大 的 原因 ， 
同时 ， 边 界 层 附 近 出 现 马 赫 反 射 现象 ， 马赫 激 波 前 
后 , 速度 由 超 音 降 为 亚 音 , 速度 波动 较 大 。 


4 Mach 
4 TKE 
/ 1200 1.5 
下 1000 忆 13 
。 800 1.1 
人 600 0.9 
400 0.7 
本 一 200 0.5 
0 0.3 
AS 01 
(a) 滑动 能 (b) 时 均 马 赫 数 云图 


图 7 干涉 区 速度 波动 


Fig. 7 Velocity fuctuation in interaction area 


图 8(a) 为 上 述 过 程 的 简化 图 ， 入 射 激 波 在 壁面 
附近 形成 马赫 反射 ， 且 因为 分 离 泡 的 存在 ， 增 大 了 
入 射 激 波 的 入 射 角度 ， 对 近 壁 区 域 激 波 反射 特性 产 
生 影 响 。 该 分 析 过 程 也 在 图 8(b) 相同 工 况 的 实验 中 
得 到 验证 。 

分 离 泡 前 ， 济 动能 较 大 ， 对 该 位 置 放大 后 得 图 
9, 左 图 显示 ,其 位 置 在 分 离 泡 前 半 部 分 右 图 为 该 
区 域 涡 量 图 , 用 流 场 中 最 大 涡 量 归 一 化 后 的 结果 , 可 
以 看 到 , 分 离 泡 的 存在 , 改变 了 其 前 半 部 分 流体 的 涡 
量 , 进而 引起 该 区 域 较 大 的 速度 波动 。 


201706.00375V1 


chinaXiv 


968 工程 热 物 理 学报 


(b) 相同 工 况 实验 结果 站 
图 8 干涉 分 析 


Fig. 8 Analysis of interaction 


(a) 干涉 简化 图 


Vorticity 
Magnitude 
0.15 


0.12 
0.09 
0.06 
0.03 
—0 


分 离 泡 出 滑动 能 
图 9 分 离 泡 处 速度 波动 针 之 


Fig. 9 Velocity fluctuation in se bubble area 


给 出 涡 量 公 式 ， 必 
FD (2) x vot (1) 
在 激 波 边界 层 作用 区 域 ， 涡 量 值 的 改变 主要 为 
斜 压 项 带 来 的 影响 分析 过 程 见 图 10。 压力 云图 中 ， 


边界 层 附近 ， 因 为 马赫 激 波 的 存在 ， 压 力 梯度 为 流 、\ 


动 方向 ， 而 密度 梯度 基本 是 垂直 壁面 方向 ， 二 
乘 不 为 零 . 做 出 斜 压 项 云图, 在 分 离 泡 前 端 ， 区 
域 , 斜 压 项 数值 很 大 ， 影响 非常 明显 , 所 以 , 分 离 泡 
前 端 斜 压 项 的 存在 , 改变 涡 量 , 引起 速度 波动 。 


Pressure | 


| | 
| J 


图 10 涡 量 公式 分 析 


Fig. 10 Analysis of vorticity formula 


关于 干涉 区 损失 , 采用 炉 生 成 率 (EGR, Entropy 
Generation Ratio) 分 析 , 公式 为 : 


1 OV 
en 二 万 
S 7 & OZ 


入 oT oT 


将 总 炉 生 成 率 分 解 为 医 性 和 不 可 道 传 热 两 部 分 ， 
公式 等 号 右 端 第 一 项 表示 由 于 入 性 带 来 的 影响 ， 右 
端 第 二 项 表示 由 于 传 热带 来 的 影响 .做 出 炳 分 布 如 
图 11 所 示 , 黏 性 和 传 热 所 引起 的 损失 分 别 用 当地 总 
炉 生 成 率 进行 单位 化 ， 表 征 医 性 和 传 热 在 炉 生 成 中 
所 占 的 比例 。 


总 依 生成 率 黏 性 部 分 
图 11 干涉 区 损失 分 布 


Fig. 11 Loss distribution in interaction area 


传 热 部 分 


激 波 作用 之 前 ， 未 发 生 分 离 ， 近 壁 区 总 烂 生成 
率 较 高 ， 边 界 导 内 除 玫 性 底层 外 ， 信 生 成 主要 由 伟 
热 不 可 逆 引 起 ; . 激 波 引发 分 离 后 ， 近 壁 区 总 烂 生成 
率 降低 ,入 插 雪 以 锋 性 为 主 ; 此 外 , 入 射 激流 和 反射 
激 波 人 山区 域 , 总 粹 生成 率 较 高 , 引起 损失 较 大 , 且 
长 要 以 竺 性 为 主 . 

-4 激 波 与 尾 迹 干 涉 机 理 分 析 

当 激 波 穿 过 尾 迹 时 ,用 密度 梯度 表征 激 波 强度 ， 
得 图 12， 和 RANS 结果 相 比 , DDES 结果 能 够 更 加 
清楚 地 反映 出 涡 结构 ,更 加 清晰 地 展示 激 波 尾 迹 的 
作用 过 程 。 

斜 激 波 穿 过 尾 迹 时 ， 产 生 反 射 现象 ， 反射 激 波 
强度 较 弱 , 很 快 耗 散 消失 ; 进一步 分 析 可 以 得 到 , 尾 
迹 还 会 改变 激 波 的 强度 ， 同 时 ， 当 激 波 进入 和 穿 过 
尾 迹 时 , 方向 都 会 发 生 明 显 改 变 . 

当 尾 迹 通过 激 波 时 , 激 波 前 后 存在 较 大 的 压力 、 
速度 脉动 ， 尾 迹 涡 结构 发 生 明显 变化 . 做 出 @ 等 什 
面 , 用 马赫 数 染 色 , 得 图 13, 尾 迹 涡 通过 激 波 时 , 整 
体 上 存在 明显 的 涡 破 碎 现 象 , 使 流 场 趋向 均匀 化 . 原 
6 的 涡 结构 , 经 过 第 一 道 激 波 后 , 已 变 为 尺度 较 小 的 
无 规则 涡 系 ,到 第 二 道 激 波 时 , 大 尺度 结构 已 基本 耗 
散 消失 。 
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图 12 瞬时 密度 梯度 


Fig. 12 Instantaneous density gradient 
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图 14 尾 迹 区 损失 分 布 


Fig. 14 Loss distribution in wake area 
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激 波 尾 迹 干涉 处 的 损失 分 析 ， 做 出 时 均 粹 生成 
率 结果 如 图 14 所 示 。 尾 迹 方面 , 尾 迹 处 粹 生成 率 较 
大 ， 但 因为 激 波 的 存在 ， 尾 迹 很 快 耗 散 ， 在 激 波 附 
近 ， 激 波 后 部 区 域 ， 因 尾 迹 带 来 的 损失 比重 已 经 很 
小 。 激 波 方面 ， 反 射 激 波 在 反射 过 程 中 存在 激 波 强 
度 地 衰减 , 经 过 尾 迹 后 , 相同 流向 位 置 时, 左 侧 激 波 
分 文 强度 方面 会 弱 于 流 道中 的 右 侧 分 支 ， 后 续 粹 增 
也 较 小 。 所 以 ， 涡 结构 的 破碎 ， 利 于 减 小 流 场 损失 ， 
激 波 强度 的 衰减 ,利用 减 小 流 场 损失 . 


3 结论 

本 文采 用 自主 开发 的 DDES 求解 器 ,解析 出 精 
细 的 非 定 常 流 场 结构 ， 为 实现 跨 音 复杂 流动 现象 的 
机 理 分 析 提 供 必 要 条 件 . 流 道 斜 激 波 在 叶片 吸力 面 
诱导 边 离 , 产生 马赫 反射 , 使 分 离 泡 前 斜 压 项 
数 y 改变 边界 层 附 近 损 失 分布 。 尾 迹 会 使 激 波 
的 反射 和 透 过 时 的 折射 , 并 改变 其 强度 ; 激 

1 后 , 大 玉 度 的 涡 结构 破碎 , 很 快 耗 散 消失 , 流 场 
趋向 均匀 化 ; 激 波 和 尾 迹 的 干涉 情况 有 利于 减 小 流 
场 损失 。 
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